'Tight Binding' methods in quantum transport through molecules and small devices: 
Prom the coherent to the decoherent description. 
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We discuss the steady-state electronic transport in solid-state and molecular devices in the 
quantum regime. The decimation technique allows a comprehensive description of the electronic 
structure. Such a method is used, in conjunction with the generalizations of Landauer's tunneling 
formalism, to describe a wide range of transport regimes. We analize mesoscopic and semiclassical 
metallic transport, the metal-insulator transition, and the resonant tunneling regime. The effects of 
decoherence on transport is discussed in terms of the D'Amato-Pastawski model. A brief presenta- 
tion of the time dependent phenomena is also included. 
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Se discute el transporte, en el estado estacionario, en dispositivos moleculares y de estado solido 
en el regimen cuantico. La tecnica de decimation permite una description completa de la estructura 
electronica. Tal metodo, en conjunto con la generalization del formalismo de Landauer, puede 
ser usado para describir un rango amplio de regimenes de transporte. Se analiza el transporte 
mesoscopico y semiclasico en el regimen metalico, la transition metal-aislante, y el regimen de tunel 
resonante. Los efectos de decoherencia en el transporte son tratados en terminos del modelo de 
D'Amato-Pastawski. Finalmente se incluye una breve presentation de fenomenos dependientes del 
tiempo. 
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I. INTRODUCTION 

With the advent of progressive miniaturization of elec- 
tronic devices, it becomes necessary to describe the trans- 
port properties of small systems within a fully quantum 
mechanical framework. As a representative example, wa 
just recall the double barrier resonant tunneling deviceu 
(DBRTD), whose current- voltage (I —V) curve presents 
a well defined peak as shown in Figjl] adapted from data 
in referencecl. This behavior, in great departure from the 
usual monotonic features shown by macroscopic samples, 
can be attributed to quantum interference phenomena, 
which arises from the coherent dynamics of the carriers. 

Let us characterize a "device", in a rough sense, as 
some specified region of material where the carriers spend 
a substantial part of the time. It might be characterized 
by a longitudinal length L x and a transverse cross section 
L y x L z = M a 2 , where a is an atomic length scale. In or- 
der to manifest its wave nature, an electronic excitation 
must propagate quantum mechanically (i.e. phase coher- 
ently) between these boundaries. Typical interference 
phenomena occurs when any of these lengths is compa- 
rable with the deBroglie wavelength X e . A condition for 
this quantum coherence is the weak coupling with the en- 
vironmental degrees of freedom within the relevant time 
scales. This condition can often be achieved by shrink- 
ing the device length scales and/or lowering the temper- 



atures until the minimum energy excitation A becomes 
large compared to the available thermal energy ksT. In 
the experiment of Fig.|l] the applied voltage is the control 
parameter for the wave length of the electrons responsi- 
ble for the transport. Therefore, the peaked I — V curve 
of the DBRTD can be viewed much as resonant peaks in 
a Fabry-Perrot interferometer. 
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FIG. 1. Scheme of a DBRTD and experimental cur- 
rent-voltage characteristic. Results adapted from Ref. 2. The 
points A, B and C in the curve corresponds to the potential 
profile shown in the inset. 

A great variety of solid-state devices developed in 
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the last decadeH satisfy the general conditions for the 
quantum manifestations described above. More recently, 
the growing interestta in understanding the conductance 
between two electrodes connected through a bridging 
molecule lies naturally in this category. In fact, dis- 
creteness of the molecular electronic energy levels is a 
clear manifestation of the quantum coherence of the elec- 
tronic states and one is left with, the question of how 
do these states transport chargeQ. The related issue 
of the propagation of a charge density excitation has 
been tackled by physical-chemists who, for a long time, 
have been, dealing with problems as diverse as conducting 
polymeraj, .charge oscillation in the intermediate valence 
compoundsQ (e.g. the Creutz-Taube ion (NH 3 ) 5 Ru-pyz- 
Ru-(NH3) 5 +5 where pyz stands for pyjjazine) and elec- 
tron transfer in photosynthetic systemaa. 

The traditional methods used to describe transport are 
not necessarily appropriate for mesoscopic and molecu- 
lar systems. The semiclassical Boltzmann equation, the 
Kubo formalism and other traditional techniques were 
specially devised to describe bulk transport, where the 
thermodynamic limit is guaranteed and the linear re- 
sponse regime is the relevant one. Hence, in order to 
deal with transport in finite size "samples" with non triv- 
ial geometrical constrains, and often in presence of non- 
perturbative fields, a new perspective had to be adopted. 
It is clear that a general Schrodinger equation must pro- 
vide the complete description. However, even in the sim- 
plest independent particle approximation, it seems diffi- 
cult to include the complex boundary conditions imposed 
by the electrodes. A way out is to treat transport of 
carriers through the "sample " as a scattering problem 
from one electrode to the other. Therefore, the trans- 
mittance T undertakes a most relevant role. This is 
essentially the conceptual framework on which the Lan- 
dauer's descriptionEl of transport is based. It involves a 
simple but conceptually new approach that inspired most 
of the advances in solid-state devices in recent years and 
is now clearing the way towards molecular electronics. 

In this article, we attempt to present a self-contained 
operational and conceptual manual based on the lectures 
given by the authors at the 2nd Workshop on Mesoscopic 
Systems. We will attempt to summarize the methods 
used to evaluate and interpret the electronic structure of 
finite systems, as well as how they are extended to de- 
scribe the coherent transport when they are connected to 
electrodes. While the effects of electron-electron inter, 
tion can become of utmost relevance in small systems 
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for didactic reasons we will not extend on this blooming 
field. Let's only mention that their understanding^ is 
built upon the methods and ideas we describe in this 
review. 

The material is presented as follows: Section II is de- 
voted to Landauer's ideas, and represents the backbone 
over which the formalism is built. Section III introduces 
the decimation method to calculate transmittances and 
electronic structure of devices and molecules. Section IV 
shows the connection between electronic-structure and 



transport obtained by Fisher and Lee. Their formula 
is derived and discussed there. Section V presents an 
overview of the known results of transport which can help 
the beginner to see the machinery in action. Section VI 
is devoted to the introduction of decoherence into the 
formalism. Section VII presents a brief outline of the 
time-dependent problem. 

A number of the most simple calculations will be 
worked out with some detail, so the unfamiliar reader 
might acquire a first feel for the topic by noting the 
simplicity of the methodology. The most advanced re- 
sults will be left merely indicated so that they can be 
reobtained by a motivated reader. Many of the exam- 
ples presented here are completely original and have not 
been published before to our knowledge. Finally, let us 
stress that we do not attempt to review the literature 
exhaustively, but simply to present our personal path- 
way through the concepts involved. In particular, the 
issues of decoherence and time dependence are the ob- 
ject of current research and the reader might still create 
his own trek. 



II. THE BASIC IDEAS OF LANDAUER 

Rolf Landauer was the first to realize that, besides the 
"sample " or device one must explicitly incorporate the 
electrodes or contacts in the transport description. He 
introduced two one-dimensional wires, left (L) and right 
(R), connecting "the sample" to electron reservoirs. The 
net current leaving the lead L is expressed in terms of the 
number of electrons available in the leads, their typical 
velocity and the probability T Rj l of transmission through 
the sample 



\L = e 



[TR.L(e)v L \N L {e)f L {e) 
-T L>R (e)v R ~N R (e)f R (e)}de. 



(1) 



The meaning of this equation is obvious: it balances 
currents. Each reservoir i emits electrons with an en- 
ergy availability controlled by a local distribution func- 
tion fi(e) = l/[exp[(e — /jJ/fegT + 1], where e is the 
energy and fa = [i + 5 fa is the chemical potential which 
is displaced from its equilibrium value p, D . The density 
of those "outgoing" states is |iV,(e) (half the total) and 
their velocity Uj. The coefficient Tr.x, computes current 
as positive provided that the particle can pass through 
the sample. Taking the i = L reservoir as reference, the 
first term is then an "out " current while the second is the 
"in" contribution. It was essential in Landauer's reason- 
ing to note that in a propagating channel the density of 
states Ni is inversely proportional to the corresponding 
group velocity: 



Ni=—. 



(2) 



2 



This fundamental fact remained unnoticed in the previ- 
ous discussions of quantum tunnelinga and it is the key to 
understand conductance quantization. Notice that there 
is no use for the traditional [1 — fj(e)] factor to exclude 
transitions to occupied final states. In a scattering for- 
mulation, the fuaal "out" states are already contained in 
the "in" statesEi Although different "in" states (e.g. 
on the left and right leads) could end in the same final 
state, unitarity of quantum mechanics assures that out- 
going contributions are orthogonal. 

To obtain the usual Landauer two probe conductance, 
one assumes time reversal for the transmittances: Tr.l = 
Tl.r. At low temperatures, the Fermi distribution func- 
tion can be safely replaced by a step function. If linear 
response can be invoked (SfXi = \i% — £f <C £f), the in- 
tegral is approximated by using the transmittance eval- 
uated at a Fermi energy. If we do not include the usual 
factor of 2 due to the spin degeneracy, one gets 

g 

I = tTr,l{Sh l - 8hr) where (S^l - 5(j, R ) = eV (3) 

from which the two-probe conductance can be calculated 

e 2 

Gr,l = -rT RiL . (4) 
h 

This conductance accounts for the "sample-lead" system. 
In a perfect 1-d conductor Tr^ = 1, and it remains fi- 
nite. In fact, resistance results h/e 2 — 25.812AT2 per 
channel, the quantum value found later in the Quantum 
Hall Effect. In contrast, the original Landauer's conduc- 
tance 

r f. P . e 2 Tr,L !k s 

G ^"i- zxl (5) 

is an attempt to extract an intrinsic property of the "sam- 
ple" . This conductance corresponds to a four probe mea- 
surement, as shown in Fig. |] (upper panel) where two 
leads are used to inject and drain the current respectively, 
while the other two additional "non-invasive" probes, de- 
noted by A and B, are used to measure the voltage drop 
in the "sample" neighborhood. They are very weakly 
coupled wires or capacitive "non-invasive" probes. The 
difference, 

R contact = 1/G^-l/G^ A ( 6 ) 

can be interpreted as a contact resistance associated to 
the lead-sample interface. 
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FIG. 2. Upper panel: Landauer's representation of a gen- 
eral electronic transport experiment. Lead at left (L) and 
right (R) inject and extract current. Weakly coupled voltage 
probes A and B at the "sample" boundaries are also repre- 
sented. Lower Panel: Portion of an actual molecular sandwich 
heterostructure between a Au-Ti left electrode and a right Au 
electrode. The bridge molecule is 4-4-thioacetilbiphenyl. The 
residue at right (X =S, Se, Te) can be substituted. 

A more general formulation of Eqs.([TJ) and (Q) for a sys- 
tem composed by many channels (e.g. spin and transver- 
sal modes) subject to different boundary conditions re- 
sults from the applicatioro of the Kirchoff law using Lan- 
dauer's conductances 

u = Ti^J [T ^ £)h{£) ~ T i,iWMte- (7) 

We do not exclude sites i = j from the sum. Here, 
the transmission coefficients may depend on the exter- 
nal parameters such as voltages, and hence accounts for 
non-linear response. The non-ohmic I — V curve of the 
DBRTD shown in Fig.0 can indeed be obtained directly 
from this formula by summing up over transversal quan- 
tum numbers (channels) in the left and right leads. Also 
A'j = t l L for every outgoing channel in the left lead and 
Mi = f-R f° r every outgoing channel in the right lead. 

Having understood the basic requirements of a trans- 
port theory one must learn to connect the transmittances 
T with the electronic properties of the "sample" and elec- 
trodes. In the treatment of solid state devices many 
researchers adopted the strategy of modeling the trans- 
port through the scattering matrices. However, when one 
deals with a molecular system, it becomes clear that spe- 
cific features of the electronic structure at the molecular 
level are relevant. Thus, at least an approximate descrip- 
tion of these properties through Hamiltonian models is 
mandatory. Let us remark that in the Landauer formal- 
ism, the "sample" is considered a finite system while the 
electrodes are considered in their thermodynamic limit 
with a continuum spectrum and acting as charge reser- 
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voirs. One then needs a formulation capable of dealing 
naturally with both situations. 



combination of atomic orbitals (LCAO) can also be ob- 
tained as a convenient approximation to the Schrodinger 
equation written in terms of the continuous coordinate x 



III. ELECTRONIC PROPERTIES THROUGH 
REAL SPACE RENORMALIZATION GROUP 
PROCEDURES. 

A. Origin of the tight-binding model 



2m 



V 2 ip(x) + U(x)ip(x) = e^(x) 



We can discrctize this equation obtaining a finite dif- 
ferences approach 



The first thing we want to point out is that tight bind- 
ing models, although naturally associated with linear 



ip(x + Ax) — 4>{x) ip{%) — ip(x — Ax) 



Ax 



Ax 



2m 



Aa 



U(x)ij){x) = eifj(x) 



(9) 



If we do the identifications 

Aa; = a: x — na\ u„ — ijj(na): 

E -U(x )■ V--£- (10) 
rj n — u \x n )y v — 2ma 2 , 

we finally obtain: 

(e - E n )u n - Vu n+ i - Fu„_i = 0. (11) 

Therefore, one is left with a discrete equation where the 
interaction is provided by the kinetic energy terms V 
which is usually short ranged. The local potential energy 
term is given by E n and, in the LCAO description, it can 
be identified with the energies of atomic orbitals. 



The above relation can be used in a "reversible" way 
since there are situations where calculations in the con- 
tinuum are simpler, as occurs in the semiclassical limit 
when the WKB approsimation substantially simplify the 
problem. 

Indeed, when the discretization is done at a scale much 
smaller than the typical atomic scale Aa; <C ao one can 
reach the molecular scale of the LCAO by a progressive 
use of the procedure that is discussed in the following 
paragraph. For the moment, we consider the tight bind- 
ing Hamiltonian: 



H = J2 E n I") H + V„, n+ i \n) (n + + V n+ i >n \n + 1) (n\ , 

n 

with (n \ip) = u n , E n = E Q and V n ,n+i = — V, from which the Schrodinger Eq.(|l^) takes the matrix form 
/ : \ (:. \f : \ 



(12) 



U n -1 
U n +1 

V : J 



E n -1 V n -\,n 
*n,n-l E n Ki,n+1 

Ki+l,n E n+ i 



U n -1 
U n 
U n +1 

V : J 



(el-H) Tt =0. 



(13) 



When the site energies are taken from a random dis- 
tribution in the range [— W/2, W/2] the Hamiltonian is 
referred to as the Anderson model and it is the standard 
model to represent disordered systems. 

In a general situation of a LCAO problem, E n cor- 
responds to any of the specific atomic orbitals used for 
each atom and the interactions parameters might not be 



restricted to first neighbors complicating the topology of 
the interaction network. 
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B. Effective Hamiltonians and Green's functions 



C. The Decimation Method 



Let's see how to construct a progressive solution of this 
problemta i^si^g the ideas of the real space renormaliza- 
tion groupt£lEZl. 

The simplest case is the two site problem 



Ul 
"2 



Ex Vi, 2 
V21 E 2 



Ul 
"2 



0, 



(14) 



whose energy spectrum is obtained from the secular equa- 
tion 



detlel-HI 



0, 



(15) 



with eigenvalues 



, Ei -+- En . 

E± = [ I 1 ± 



Ei — E 2 



V h2 V 2A . (16) 



An alternative procedure is to write the linear equation 
explicitly 



a) E1U1 + Vi. 2 u 2 = eui, 

b) V24M1 + E 2 u 2 = eu 2 . 

From b) u 2 = V 2 ,\ e ^ E Ui, and substituting in a): 



(17) 



The preceding paragraphs have introduced the basic 
tools by noting that the presence of other sites has the 
effect of shifting or splitting the energies of the site we are 
looking at. Let us see another example, a Hamiltonian 
with three sites: 



El Vl,2 Vl,3 

V 2i i E 2 V 2 ,3 

VJj,! 1/3,2 E3 





j Ul 


u 2 


= e u 2 


U3 ) 


\ U3: 



(22) 



where we have included second neighbor overlaps. The 
idea is to reduce this set of three linear equations to a 
smaller problem as shown in Fig.^. There are many sit- 
uations where the physics and chemistry indicates that 
this is a physically meaningful model. We obtain the 
value of u 2 from the second row above 



"2 



:-E; 



V-2 



■ui 



2.:i 



E 



■u 3 , 



(23) 



and introduce it into the first and third row to obtain 
two coupled non- linear equations: 



Ei + Vi 2 - 



1 



V 2 ,i ui = eui. 



e-E 2 

The eigenvalue is obtained from the condition 

1 



(18) 



e - [Ei+ V 1<2 



e-E 



■V2,l\ 







(19) 



Ai(e) 



The reader can check that both exact eigenenergies, can 
be obtained from this equation. The second term in the 
parenthesis has a clear physical meaning if it is identi- 



fied with an "effective potential" Ai(e) 



\Vl,2\ 2 

s-E 2 



which 



corrects the non-interacting energy Ei of the site as 

Ei=Ei+Ai(e). (20) 

For an asymmetric case, \Ei — E 2 \ > \Vi y2 \ , one can ob- 
tain a good approximation to the corrected site energies, 
by doing the evaluation at the old eigenvalue: 



Ei ~ Ei + Ai(£/i), 



(21) 



which becomes equivalent to the second order Rayleigh- 
Schrodinger perturbation theory. That is, as shown in 
Fig.|3|, we have gone from the two orbital problem to a 
single effective orbital with a "dressed energy " . 



Ei V1.3 m 
^3,1 E 3 )\u 3 



ui 
U3 



where 



E 1 (e) = E 1 +V 1 , 2 - 



1 



-V 2 ,i, 



' e - E 2 

i?3(e) = E3 + V32 ^-^2,3, 

e — E 2 

Vl,3 - V Xt3 + Vi, 2 — V2.3, 

S — £j2 



and the notation is self-explanatory. 



Vi, : 



(24) 



(25) 



FIG. 3. Construction of the self-energy correction. 
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Ei E 2 E 3 

V 1/2 V2,3 



V 2 ,i 



Ei 



V 



3,2 



Vi, 3 



V 3/] 



D. Effective Hamiltonians and Green's functions 



The Green's function provides an alternative frame- 
work to discuss the solutions of the Schrodinger equa- 
tion. Besides the intuitive structure of thejierturbative 
calculations in terms of Feynman diagrams^!, they have 
the additional advantage of a clear connection to trans- 
port properties, and finally a special role in the Quantum 
Field Theory which allows a definitive treatment of the 
many-body problems. 

The retarded (advanced) Green's functions matrix is 
defined as a function of the complex variable (e ± i?7) in 
the upper (lower) complex plane. 



G R{A) (e) = [(e±i7?)I-H]- 
Let's evaluate a diagonal element: 



(29) 



FIG. 4. Construction of an effective Hamiltonian through 
decimation. 

Therefore, the procedure of elimination of variables is 
very general and can exhaust the degrees of freedom of 
the finite system, providing a systematic way to reduce 
the dimension of the Hamiltonian. This is done at the 
cost of transforming the linear equation into a non- linear 
one which, however, can often be linearized in the region 
of interest as in Eq. (pjl|). The general recipe valid for a 
string of interacting atoms requires the calculation of: 



• Ai( n )j the energy correction to the 1 st atom when 
sites until atom n, to the right (+) inclusive, are 
eliminated. 

• A~ +1 (e), the energy correction to atom (n + 1) 
when all the atoms at the left (-), except the 1 st , 
have been eliminated. 

• Vi.n+i, the effective interaction among layers. 



We resort to the recursion formulas: 

f 



^5„)(s) = A+ (n _ 1) +V 1 , n 



£ - E„ - An 



14,i 



A B+1 (e) = K+i,-. 



f 



— V n ,n+1 



e-E n - A;. 

^l,n+l = V\,n — -J—:Vn,n+l + ^l.ra+l 

£ — E n — A n 



(26) 
(27) 
(28) 



Notice that we have not restricted the iterations to near- 
est neighbors. From the point of view of the perturbation 
theory, it can be checked that the decimation procedure is 
equivalent to an exact summation of all the perturbation 
orders of the Wigner-Brillouin series and hence equiva- 
lent to the inclusion of all Feynman paths that start at 
layer and end at layer N + 1 . This will be best viewed 
in the language of the Green's function described in the 
next subsection. 



s + irj-Ex- V-i 



1 



1,2 



e + irj - E 2 



Vl, 



iTi 



1 



e + IT] - Ei 

Therefore the renormalized site energy is as before 

1 



(30) 



Ei — E\+ Vi t 2 



e + ir) — E- 



(31) 



]f = A x - iTi 



The local densities of states represents the weight of the 
exact eigenenergies on the old states, which we can check 
for our simple two state model, is: 

JVi(e) = \u h+ \ 2 S(e - E + ) + \u h _\ 2 6(e - E_) (32) 

,2 1 V 



lim <j | mi 



7T (e - E+f + Tf 
2 1 V 



it (e-E-) 2 +T] 2 J ' 
where 6(e — E±) are the Dirac delta functions at E± 



(33) 



$[(E 1 +E 2 )±hu], with hu> = y/(E 1 ~E 2 ) 2 +4|V r i, 2 | 2 . 
The eigenstates are |+) = |1) + u 2 ,+ |2) and |— ) = 
iti._ |1) + «2.- |2) with the coefficients 



Vi 



Ul,± 



1.2 



2\Vi, 



1± 



Ei — Eq. 



1/2 



"2,± 



± 



Vi, 



21V, 



1.2 



it 



Tulo 



Ei — E 2 



1/2 



(34) 



(35) 



G 



An immediate advantage of the Green's function formal- 
ism is that although it is obtained through a finite num- 
ber of algebraic operations (matrix inversion) it contains 
information on both the eigenenergies and eigenfunctions 
which involve the transcendental operation of finding the 
roots of a polynomial. In particular, any eigenvector 
component is obtained from the generalization of Eq. 
(E32I) which results: 



Ni(s) = -- lim ImGfAe + irj) 

7T »?^0 



(36) 



The connection with the infinite order perturbation 
theory is immediate just expressing the denominator in 
Eq. (|30j) through its series expansion 



irj-Ei- Vi 



1 



1.2 



s + irj- E : 



■V li2 



Vr,2G^l R {e)V 2 , x 



_ n {o)R 



n=0 



Vi, 2 G$ R (e)V 2A 



(37) 



(38) 



(39) 



- G[T(e) 

+G{i%)V 1 , 2 Gii%)V 2 , 1 G{°l\) 



= G$ r (e) + G { °{ r (e)X 1 , 1 (e)G* (e) 



(40) 
(41) 



The last line corresponds to the Dyson equation while 
the series expansion is represented diagrammatically in 
Fig.|. 



/\ 



+ 



* + 



A 



" + 



FIG. 5. Feynman diagrams in tight-binding representation. 
Double and simple lines represent the exact and unperturbed 
Green's function respectively. Hopping interactions are rep- 
resented by the dashed lines that modifies the indices of a 
Green's function. Two conjugate interactions are associated 
by the dot into a single rung. 

In systems with a finite number, N, of states Gjj-(e) 
is a well behaved meromorphic function except at the N 
poles. In this condition, it is satisfied that 



finite system 



T(e) = -Im[E( £ )] = 0, 
E(e) = A(e) = Re[S(e)]. 



(42) 



Therefore, by considering real energies we can drop the 
supra-index that distinguish the retarded and advanced 
Greens function. 

Let us come back to the three site model. We may sum 
up the infinite perturbation series to obtain 



Gi,i(e) 



(43) 



e-E 1 - Vl, 



-V 2 a - Vi, 



e - E 2 + V 2 



2.3 



E-E a 



-V, 



e-E- 



■V 3 ,i 



:',:2 



or equivalently use Eq. ([25|) to regrouping the terms 
into: 

1 



GiAs) 



E1-V1 



1 



(44) 



1,3" 



En 



-V 3 A 



which gives the obvious connection with the real space 
decimation procedure. 

Many other useful recursion formulas can be obtained 
by similar procedures. For example, in a system with 
states {0, 1, 2, 3, • • -N, N + l}, the effective hopping is re- 
lated to the non-diagonal Green's function in the isolated 
bridge with states {1,2,3,- • • N}, 

-1 



N+l 



V . 



N 



E 



N 



Vn,n+i, 



(45) 



— Vo,i yGi t ff—iVff—i t ff [e — En — Ajy] jv^jv+i, 
= V 0ll {G liN }V N , N+ i. (46) 

This equation is complemented with the ones for the self- 
energies from sites at right 



Ao(jV)( £ ) = ^0,1^1.1^1.0) 

and at left 

^jv+i(jv)( £ ) = Vn,n+iGn,nVn,n + i. 



(47) 



(48) 



One can then calculate by iteration of the non- 
diagonal Greens Function in progressively bigger systems 
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{0, • • • N + 2}, {0, • • • N + 3}, and so on, with an equation 
that is obtained equating the terms between brackets in 
Eq.©: 



,(JV+1) _ n {N) 
r l,AT+l — u l,jV 

which is used in conjunction with 



f-y(l\)-rr /~<{N + 1) 

U 1,N+1 — U 1,N V N,N+1^ N+1 , N+1 , 



(49) 



G 



(JV+l) 
JV+l.JV+l 



el — En + i — Vn+i,nG^nVn+i,n 



(50) 



the continued fraction expansion of the diagonal Greens 
function. 



E. Going beyond one-dimensional systems. 

The decimation scheme can be generalized to any di- 
mension as long as we proceed in a "layer by layer" elim- 
ination. Every site in the above procedure now becomes 
an n x n-matrix where n is the size of the layer. In Fig. 
^ we consider a square lattice of orbitals forming a strip. 
Remembering that Im[S] = we write the general equa- 
tion for the real self-energies. We proceed to eliminate 
one by one the layers from 1 to the (N — 1). Adopting 
the matrix notation A^( n ) (e) is the self-energy correction 

to the 1 st layer when all layers (to the right) including 
the n th have been eliminated. A~ +1 (e) is the self-energy 
correction to layer (n + 1) when layers to the left have 
been eliminated. Vi.„ + i is the effective interaction be- 
tween layers. These are again evaluated with the iterative 
procedure 

Je) = A+ + Vi, n — — ^ — V nil) (51) 



E, 



A- +1 (e)=V„ 



+ 1,71 



el - E„ 



:V„. 



71+1, 



Vi 



71 + 1 



Vi, 



1 



£l - E„ - A n 

with the initial values: 

A+ (1) (e) = 
A 2 -( £ ) = 0, 

V 1 ,2=V 1 , 2 . 



- V„,„ +1 + Vi, 



71+1, 



(52) 



(53) 



(54) 



The expression one has to evaluate are expressed as Ma- 
trix Gantinued Fractions, which are numerically very 
stableEJ. Notice that the elimination of intermediate lay- 
ers produces effective interaction among the sites in the 
first layer, i.e. it modifies the intra-layer interactions as 
prescribed by the non-diagonal elements of A. Inter- 
layer interactions are always contained in V. Again, the 
last term in Eq. (|53|) allows the consideration of interac- 
tions going beyond nearest neighbor layers. 




n+l n+2 N 



.0— f f" 



— o 

*"•<# — i 



FIG. 6. Scheme of the decimation of a finite tight-binding 
strip. As the elimination of intermediate layers proceeds the 
effective site energies and interaction appears. 



F. Decimation in Molecules 



It is clear that the described procedure i.s .very well 
suited for application to molecular systemst3c3. Just to 
fix ideas consider the p z orbitals in the carbon backbone 
of an organic molecule represented in Fig.j^. Assume that 
we are interested in the study of how charge can be trans- 
ferred from site L to site i?, a typical problem in pho- 
tosynthesis and molecular electronics. Circles represent 
the 7T orbitals with given local energies, lines are hopping 
interactions which produce the electron derealization. 
One might start decimating the dangling ends. Sites en- 
ergies in the back-bone are then renormalized. Next, we 
identify the branching nodes and eliminate the bridging 
sites, this gives new site energies and an effective hopping. 
Finally one eliminates all the remaining bridging struc- 
ture obtaining an effective two orbital molecule. Since 
the procedure is exact one obtains the exact spectrum 
independently of the election of these orbitals. However, 
in transport one wants to keep the atoms that have ma- 
trix elements that enable the transfer of charge with the 
"external world" . This will become clearer later. 
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CO 



FIG. 7. Sequence of the decimation of a molecule into an 
effective two site problem. 



d) The resulting highly singular resolvent of a big sys- 
tem through matrix inversion is in general numerically 
unstable. In contrast, the decimation, is a numerically 
stable procedure. The natural instabilities of the system 
such as exponentially growing eigensolutions are explic- 
itly used in favor of the convergence of .the method. For 
further discussion on this issue see Ref.t^l. On this basis 
one establishes that Eq. (|52| ) which defines a Matrix Con- 
tinued Fraction procedure is more accurate and stable 
than Eq.(|50|) which uses a different iteration procedure 
to calculate an equivalent magnitude. 



The spectral and transport properties of model 
molecules are further discussed in Ref.flla). Let us men- 
tion that the procedure allows to visualize a situation 
of maximum coupling when El = E^ which is a res- 
onant situation. This paper also discusses for the first 
time a situation of maximum decoupling which occurs 
when Vl,r — and was named "antiresonance" or mini- 
mal effective decoupling of the centers. This phenomenon 
is caused by the interference between different pathways 
when energies of the pathway molecules lie between en- 
ergies of the other. This generalizes, to the transport 
case, a concept introduced by FanoEil in ionization spec- 
troscopy. 



G. Advantages of the decimation procedure 

One might think that instead of decimating a particu- 
lar structure one might invert (el — H) directly. However, 
there are various advantages in favor of the decimation 
procedure explained above: 

a) The scheme adapts naturally to increase the size 
of the system without having to recalculate the matrices 
again. 

b) The sparse nature of the Hamiltonian is naturally 
included in the recursion formula saving storage memory 
and iteration steps. 

c) The decimation scheme has implicit a deep physical 
insight of the system and reveals the self-similarities of 
the systems whenever those properties are present. 



H. The Thermodynamic limit: Dyson Equation 

We have seen that any finite Hamiltonian matrix of 
finite size can be solved through the decimation method 
obtaining a set of discrete eigenstates. One might won- 
der what happens when the number of orbital is actu- 
ally infinite. How and when does the continuum spec- 
trum appear?. Let's come back to the one-dimensional 
model. The decimation procedure allows to deal with a 
simple but already non-trivial case: that of an ordered 
semi-infinite chain. In this situation: V^n+i = V and 
E n = Eq for every n. The Dyson equation is: 



£ n = V- 



1 



E 



-V. 



-*n+l 



(55) 



and must include the simple fact that every site sees to 
the right an infinite chain: £„ = X for every n. That is 
a new way to present the Bloch theorem. 



V- 



1 



E 



-V = A=Fir. 



(56) 



The striking fact is that even when we are working with 
real quantities,-the solution of this equation may lay in the 
complex plane£3. There are two possibilities for the imag- 
inary part. We call retarded self-energy to that which 
would cause a decay in the time evolution of the wave 
function amplitude. The solution of the second order 
equation gives: 

J 



A = 



E-Eg 

2 



e-Eg 
2 



2 

£~Eg 



' -V 2 i£e-Eo>2\V\, 
if \e-Eo\ <2\V\, 



(57) 



and 



r = 




o 



e-E 



if \e-E \ >2\V\ 
if \e-E \ <2\V\ 



(58) 
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We see that in the region were the Block solutions oc- 
cur, i.e. when e = Ek = E — 2Vcos(ka), is the re- 
gion where the spectrum is absolutely continuous, one 
has r(e) 7^ 0. In the region of the real axis, the Green's 
functions indetermination represents a branch cut, and 
it becomes non- analytic at the spectral support. How- 
ever, the infinite dimensionality of the Hilbert space is a 
necessary but not sufficient condition for the continuum 
spectrum. For example, for 1-d systems with disorder in 
the site energies described by the Anderson model, the 
self-energies are always real even when the system is infi- 
nite i.e. there is genuine "phase transition" in the nature 
of the states when one goes from finite systems to infi- 
nite ones. In an infinite system, disorder can produce 
a localized-extended transitioned first described by P. 
W. Anderson and commonly referred as metal-insulator 
transition (MIT). The thermodynamic limit that makes 
possible the study of such transition is to study the ob- 
servable O (decay rate, density of states, etc.) in the 
limit 



lim lim O. 

r)— >0 N—*oo 



(59) 



The order of this limit was implicit when we searched for 
complex solutions for the Dyson equation by imposing 
the "Bloch theorem" . 

We must remember that the MIT is not exclusive to 
the disordered systems. It is also of frequent occurrence 
in ordered-systems under the action of incommensurate 
potentialstd For example a 1-d systems where the site 
energy is of the form E n = W cos [Qua] with 2n/Q in- 
commensurate with the lattice length a provided that 
W >W C = 2V. 



I. Density of States in unbounded Systems: Crystals 

In the case where the chain extends to both sides of 
site one gets two contributions to the self energy cor- 
rection. All diagonal terms of the Green's function are 
identical 



1 



e-So-2S- r ' 



e — -fro — I 2 



(60) 



from which the density of states per site (DoS) can be 
evaluated, 



N [d =i] (e) 



1 



1 



2ttV 



(61) 



2V 



which has the characteristic van Hove singularities of one 
dimensional systems at the band edges. 

The fact that the imaginary part in the self-energy is 
associated to the "irreversible decay" of the state can 



be formalized by studying the Fourier transform of the 
Greens function to obtain the return probability 



P , Q (t) = |(0|exp [-iHt/h]\0)\* 

Gq (e) exp [—iet/h] de 



\Mtv/h)\ 2 ~ ~(tv/h)-\ 

7T 



(62) 
(63) 



where Jo is the Bessel function. For long times, it repre- 
sents an irreversible "super-diffusion" , and hence differs 
substantially from a classical random walk process in an 
infinite chain. 

Densities of states in higher dimensions for hyper-cubic 
lattices can be obtained from those in one dimension by 
ussing the fact that variables are separable. This allows 
us to split the energy contributions of each dimension 
E(k) = E + E(k x ) + E(k y ) + E(k z ), leading to a con- 
volution expression for the DoS in dimension d=2 and 
d=3: 



(64) 



N [d+1] (e) = J N [d] (e - e>) N {1] (e 1 ) de'. (65) 
For example, in a square lattice we get 

(66) 



and in general 



1 



N m (s) = -K 



1- ( — ) 2 

y 2V' 



where K is the complete elliptic integral of the first kind. 
Fig._| show the qualitative features of these DoS. See 
Ref.Effor a more extended discussions on these points 
and plots of the DoS in higher dimensions. 





-2V o 2V £ -E» -4V 4V 

N(£), 



d=3 



-6V 



e-E„ 



FIG. 8. Schematic Density of States in 1-d, 2-d and 3-d 
systems. 
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Now that we know how to evaluate some basic self- 
energies, we can calculate the Green's function in a great 
variety of model systems. 



J. Surfaces in a semi-infinite chain. 



We leave to the reader the evaluation of the local DoS 
both numerically and analytically, as well as the calcu- 
lation of the condition for the appearance of localized 
states and resonances. For V s = V, a localized state may 
appear. It shows up as an isolated pole in Gf; s (e) in the 
real axis at 



Let us consider the semi-infinite chain {s, 1, 2, 3, • ■ ■}. 
The Green's function at the surface site is 



"•1,1 



1 



(67) 



where we have used the subindex 1 for the Green's func- 
tion to stress on the fact that we count orbitals starting 
at the surface. One gets the density of states for the 
surface of the chain 



1- 



e-E 
2V 



(68) 



Incidentally let us note that an identical DoSris obtained 
when the Hamiltonian is a Random Matrixes. In this 
case one can use the Lanczos method to tridiagonalizc 
the matrix in the infinite dimension limit and see that it 
corresponds to the ordered chain we just discussed . 

The return probability to the surface site, as defined 
in Eq.(|62j) gives 



^1,1 (*) = 



—j x {2Vt/n) 



\ 1 


' n ' 


7T 





1 3 



(69) 



indicating a decay of the surface state. This exact result, 
which differs from the usual exponential and diffusive 
laws usually adopted for the empirical description of de- 
cay phenomena, describes an "irreversible" escape from 
the site. This contrasts with finite-systems where recur- 
rences, called mesoscopic echoesc3E3, appear at a typical 
time estimated as h/A, called the Heisenberg time. 

Perhaps, the most important conclusion to be drawn 
from the previous results is the correction to energies of 
finite systems to account for its contact to infinite sys- 
tems. 



E s = [E. 



E Q f + V 2 



I [E s - E Q 



(71) 



provided that \—Eq\ > 0. 

One can also check that for \E. 



E„\ < V and 



< 



^0 1 ^ y | y 

1, a resonant state will appear a pole in the complex 
plane with a real part overlaping with continuous band. 
Its energy is: 



E, 



E s + 



A(E S 



(72) 



The return probability is well approximated by an expo- 
nential law as prescribed by the Fermi Golden Rule: 



P s , s (t) ~ exp[-t/r s ], 



(73) 



with l/ Ts = ^\V s \ 2 N sulf .(E s ) ~2T s (E s )/h. (74) 

Again, one sees that the time dependencies that appear 
in Quantum Mechanics can be very different from those 
present in a classical random walk. 



L. Localized state in a branched circuit. 

There is a highly non-trivial result@ that can be easily 
obtained with the help of the above formalism. This is 
the existence of localized states in branching regions. 

In general, for a site x to which z semi-infinite chains 
are connected, the Green's function at the crossing site 
between branches can be written as 



G 



R 

x.x 



En 



(e - E ) - iJW* - (e - EoY 



K. Adatoms and Surface states. 



(75) 



The method given above allows us to extend the deci- 
mation procedure for finite systems described in the pre- 
vious section to composite systems (finite +infinite) . Let 
us discuss the simplest but still non-trivial example. In- 
deed one can think of the 0-th site as an adatom in the 
surface of a metal. In that case its energy E s and hopping 
element V s are different from those of the bulk, Eq and 
V respectively. The Green's function evaluation leads to 



U s,s — 



1 



Ea — U* 



{e-E )-\^/lV*-(e-E ) 2 



(70) 



Which besides of the expected branch cut in |e 
2V , it presents isolated poles at 



E 



branch 



= E ±2V 



4(z-l) 



E \ < 



(76) 



These solutions, are not only present in a branching poly- 
mer, but in a conveniently grown heterostructurecJ. The 
existence of a topologically confined states in the casa 
z = 3 hap-become the basis for quantum transistorsEEl 
and lasers^. 



11 



M. Representation of the environment through 
self-energies 

First we recall that one, can. always eliminate micro- 
scopic degrees of freedomtSiIa generating an effective 
Hamiltonian that accounts for them exactly. This pro- 
duces effective interactions and energy renormalization 
which depend themselves on the observed energy. Fur- 
thermore, one can include the effects of a whole lead in 
a Hamiltonian description through an imaginary correc- 
tion to the eigenenergies. In fact, an electron originally 
localized in the region called "the sample" should even- 
tually escape or decay toward the lead. Hence, there is 
a escape velocity associated with the energy uncertainty 
of a local state i: 



2a a 
"ST » = ' 



(77) 



Other interactions that one usually considers produce de- 
cay rates reasonably well described by the Fermi Golden 
Rule: electrons in an excited atom decaying into the 
continuum, or propagating electrons decay into differ- 
ent momentum states by collision with impurities pro- 
ducing phonons or photons. In some of these cases we 
have to add some degrees of freedom to the sum (the 
phonon or photon coordinates). A process a may pro- 
duce a complex contribution a Yj R which adds into the 
total self-energy Y, R . In particular this is the interpre- 
tation one must give to the small imaginary part 77 we 
introduced in the definition of the Green's function of 
Eq. (p9|). The effective Hamiltonian becomes: 



interactions 



with 



S«W=^(A i Tir i )(H 



(78) 
(79) 



That is, we can not find cigenstates because the complete 
Hamiltonian is not separable into a product of states of 
the unperturbed system and those of the environment. 

In terms of the Green's function the resulting equation 
is, in the matrix representation: 



g r = G (o)R , Q(o)i? s fiQfl 



(80) 



The diagrammatic expansion for the perturbed Green's 
function is shown in Fig. |^. 

Although the Green's function formalism seems to in- 
troduce some extra notation, this has various conceptual 
advantages. For example, it is straightforward to use 
Eq. ( |sb"| ) to prove a version of the optical theorem which 
relates the imaginary part of the forward scattering am- 
plitude with the total cross section of the pertubation. 
One starts by writting Eq. ( p0| ) in the form: 



( G (°)«)-i _ (G R )-\ 



(81) 



and an equivalent one for Y, A . The difference between 
them is: 



- S A = (G A ) 
= (G R ) 

from which we obtain: 



-1 [ G * - G A ] (G A ) 



[G R - G A ] = G R [X R - S A ] G A , 



(82) 



(83) 



which has deep physical significance. It is an inte- 
gral equation relating the local densities of states given 
Eq.(^6) and the decay rates provided by Eq. (|77|). Per- 
haps, the most important advantage of Green's func- 
tions, is-that they can be used also in the Quantum Field 
TheoryExEJ to deal with the many-body case and are a 
relevant tool in the clarification of the problem of tun- 
neling time. 



/ » 
- + - 



k I 




phonon 
emission 



exploration 
to the leads 



FIG. 9. Diagrams shown before can be rearranged into 
the self-energy correction E. This equation is valid for in- 
dependently of the approximation used to calculate E. In 
the lower panel self-energies obtained through escape to the 
leads and electron-phonon interactions are shown. The es- 
cape self-energies contain the hopping (dot) and a propaga- 
tor in the lead. An electron-phonon self-energy is evaluated 
in terms of the sample's electron (straight-line) and phonon 
(wavy-curve) Green's functions. 



IV. THE FISHER AND LEE FORMULA 



)iH 



The original demonstrationcj of the connection be- 
tween transmittance and Green's function is based on 
properties of the scattering matrix in the continuum 
model. It is also presented in a simplified form in the 
book by DattaEJ. Here we are going to present the main 
lines of a demonstration based on properties of a molec- 
ular orbital model in its decimated tight binding version. 

The isolated sample can be decimated into two sites 
getting an exact non-linear effective Hamiltonian: 
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Ho = E 1 11) (1| + E N \N) (N\ + V 1>N |1) (N\ + V N>1 \N) (1| 



(84) 



where the parameters E and V are functions of the variable e. 

We attach electrodes or leads at each site. Both leads are perfect semi-infinite 1-d chains where the dispersion 
relations are 



E k L = E L + 2V L cos[A;lo], 



from which the group velocity is 



vl = 



1 dE, 



kh 



2aVL ■ \U 1 



Analogous relations for E^r and v R hold in the right channel. When the perturbation 

Wo-icads = V tL (\l) (0| + |0) (1|) + V tR (\N) (N + 1| + \N + 1) <JV|) 
is turned on, the new single particle eigenfunctions components 



Uk(an) = < 



' exp[ifcna] + r exp[— ikna] 
A 
B 

texp{ikna] 



if n < 0, (na e L), 
if n = 1, 
if n = N, 
if n > iV + 1, (na e R), 



(85) 



(86) 



(87) 



(88) 



are determined from the stationary Schrodinger equation i^el — T-Lj \ip) = 0. We have only four non trivial terms 
which contain either A or B. One eliminates A, B and t and evaluates 



Tr.l = l- r*r 



= 4 



V tL 



Vr 



d^lsm^a])!^ 



N 



V t R 



Vr 



(\V R \sm[kRa}) /[dd*] 



(89) 



(90) 



where 



d 



e-E 1 + 
x e - E N + 

- Wi, N v N ,i] . 



V tL 



v L 

V t R 



Vr 



{\Vl\ cos[k L a] - i \Vl\ sin[fcLa]} 



{\V R \ cos[k R a] - i \V R \ sm[k R a}} 



In order to compare this with the Green's function we 
remember that the effective site energies include the self- 
energies produced by the leads 



and 



v t 



LL 



V L 



(\V L \cos[k L a\ - i\V L \sm[k L a\) (91) 



= L Ai -i L r! 



V 



tR 



Vr 



\V R \cos[k R a] - i|VR|sin[fc fl a]) (92) 



r a 



N 



R r 
1 N . 



If the lead variables are decimated they produce an ef- 
fective potential of the form 

Weff. = Ho + L Si |1) (1| + R Z N \N) (N\ . 

This maintains the structure of an effective two site 
Hamiltonian. From this it is simple to compute the four 
components of the exact Green's function of the system 



G?„ - <1| 



el-H 



off. 



I TV) 



(93) 



V hN 



[s - (E 1 + L Ef)] [e - (E N +« E*)] - V hN V NA 
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obtaining the Fisher-Lee formula 



T R , L (e)=r-) v R v L Gl L (s)Gt R (s). (94) 

Originally, Fisher and Lee considered only the escape ve- 
locity to the leads (i.e. a = (3 = lead). D'Amato and 
PastawskEj were the first to realize that because of Eq. 
( |S6| ) one could write transmittances in terms of the imag- 
inary part of the self energies. Immediately, our point is 
that any other process which contributes to decay giving 
an imaginary contribution to the self-energy would be 
described by Eq. (|94|). In particular, this will be true 
for a "decoherence" velocity as would be the case of the 
electron-phonon rate described in the previous section. 
In a modern notation: 



channels connected to their respective reservoirs. When 
the two wires L and R have finite cross section, each 
wire contains an integer number, L M^ £F - ) and fl M( £F ), of 
propagating (one-dimensional) channels (i.e. transversal 
modes) allowed by the Fermi energy. Each electronic ex- 
citation leaving the wire L in the propagating channel i 
with velocity vn ^ 0, enters the wire R in the propagat- 
ing channel j with velocity v R j ^ 0. The conductance in 
the linear response regime is 



(97) 



2 M "M 



This can be written by combining Eqs. ( |95| ) and (|97|)into 
a compact matrix notation: 



T aR Me) = [2 a r R (e)} G^ L (e) [2^T L (e)] G$ L>aR (e). Gr l = ^ Tr [T R {e) G% L {e)T L {s) Gl R {s)} 



(95) 

The left supra-index in T indicates the process or channel 
associated to the electron decay from the spatial region 
identified with the right subindex. Both indices appear as 
sub-indices in the Green's function. Notice that °Tl(£), 
which corresponds to the interaction with the "source" of 
particles, was arranged between two Green's functions. 
a T R (s), representing the properties of the "sink", was 
placed at left. This order is important in the matrix rep- 
resentation. 



A. An example: branched circuit. 

In order to apply the Fisher-Lee formula to a sim- 
ple but highly non-trivial example, let us consider again 
the circuit with z equivalent branches. One might think 
that a plane wave incoming at the node from one of the 
branches will have a probability \j(z — 1) to be trans- 
mitted to each of the others, i.e. that it will behave as a 
perfect (z — l)-splitter. The true answer comes from the 
quantum evaluation giving 



T = 



1 



l + (^) 



(e-Eo 



(96) 



{AVf - (s - E a f 



T has a maximum of 4/z 2 , lower than the classical value. 
Additionally, the transmittance goes to zero as the energy 
approaches the band edge. This is due to the presence of 
the localized state found in Eq. ( |76| ) which precludes the 
occupation of the branching region by the propagating 
states. 



B. An extension: the multichannel case. 

In fact, one can see that it is not difficult to gener- 
alize Eq. (0) to any number of incoming and outgoing 



h 

The sum of initial states at left is the result of the prod- 
uct of the diagonal form of the broadening matrix r^(e), 
while the trace corresponds to the sum over final states. 



V. APPLICATIONS OF LANDAUER'S 
COHERENT CONDUCTANCE 

A. Lyapunov Exponents and Localization Length. 

The calculation of actual eigenstates in disordered po- 
tentials of an infinite system is a very subtle mathemat- 
ical problem. In contrast, the computation of Green's 
functions and conductances of arbitrary size can be done 
at a relatively low computational cost. The sparse nature 
of the single particle Hamiltonian when it is expressed in 
a localized basis (tight-binding approximation) allows the 
implementation of a recursive calculation of the Green's 
function. These methods, originally devised to deal with 
inhomogeneous and disordered systems, allow the cal- 
culation of the quantum transmittances such as T R l. 
Typically, it requires L x /a operations (products and in- 
versions) with matrices of size M x M. 

Historically, Landauer's reasonings were decisive in the 
development of thc—thriory of electronic localization in 
disordered systemsEflO. In particular, the method de- 
veloped above can be used to proof rigorously that in 
one dimension there are no extended states. Every eigen- 
state should decay exponentially away of some particular 
localization center. The exponential rate corresponds to 
the Lyapunov exponent 7 of the iterative Dyson Equa- 
tion. The associated length, £ — 1/7, is called localiza- 
tion length. 

The Lyapunov exponent can be evaluated from the 
Green's functions G\ iTL and G n . n using the limit relations 



7(e) = — lim 
a L— >oo 



Gix(e) 



Gi,l(e) 



(99) 
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This definition is completely general and works either 
for ordered and disordered situations. To check how it 
works, let us first consider the 1-d ordered chain. The 
solution of Eq. (|99| ) is analytic 



7(e) 



-arccosh 

a 

o 



E 



2V 



for \e-Eo\ > \2V\ 
for \e-Eo\ < \2V\ . 



(100) 



This clarifies that "/(e) is just the analytical continua- 
tion of the wave vector k(e). Already in this case, we can 
identify two regions: one which coincides with the Bloch 
band where localization length is infinite and the rest, 
where the states, if any, are localized. This is the region 
where the described states in the branching region and 
the surface states lie. 



B. Localization in a Conjugated Polymer. 

In order to consider a more realistic example of specific, 
interest in molecular electronics, let us consider a modelEI 
for the simplest conjugated polymer: polyacetilene. This 
can be described by the Hamiltonian, 

Ho= {^a|n}<n + l| + Vi|» + l}<n + 2|+c.c.}, 

n(odd) 

(101) 

where V 2 = — (K- + 5) and V\ = —(V n — 8) where V* > 
is the energy associated with a tt — it bonding. This is 
represented schematically in Fig.|l^. S is the additional 
energy involved in the dimerization to form the alternate 
double bond (conjugated system). Hence, the unit cell 
has two C-H monomers. Again, the Lyapunov exponent 
can be calculated either analytically or numerically giv- 
ing 



- arccosh 



'-(V, 2 +v 2 2 ) 



2KV, 



at the gaps 



7(e) = \ Oat the valence band, V 2 + Vy < e < V 2 - Vi 

Oat the conduction band, — V2 + V\ < £ < — V<i — V\. 



(102) 
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FIG. 10. In the upper panel the polyacetilene macro- 
molecule is represented. Double bonds and simple bonds al- 
ternate. This gives an alternation in the energies associated 
to 7T-7T bonding. The lower scheme is the Lyapunov exponent 
as a function of the energy. Scales are discussed in the text. 

This simple situation hints at the solution of a the more 
complex disordered systems. Conceptually, a disordered 
sample can be visualized as a periodic array of disordered 



unit cells of size L = Na, with a progressively increas- 
ing N. Then, the localization transition appears as the 
dominance of the "gaps" over the "bands" whose width 
(spectral support) decreases as V exp[— 7JV0]. With the 
decimation methodology developed above, the numerical 
evaluation of the Lyapunov exponent is straightforward. 

C. Metal-Insulator Transition. 

In one dimensional systems, the exponential decrease 
of the conductance follows directly from the statis- 
tical properties of the transmission probability T ~ 
cxp [— L x /£\ in a disordered system. Therefore, the 'four- 
probe ' Landauer's conductance given by Eq. (||), by de- 
scribing the intrinsic properties of the "sample" , contains 
the correct scaling behavior of the conductance ranging 
from the non-extensive behavior when L x > ^, 

Gfl P L(^)^^exp[-L x /e], (103) 

to the expected Ohmic behavior described by Drude's 
law in terms of the mean free path I = £/2, 

Gfl P i(ix) = e 2 N v F £/L x , (104) 

when L x < I. Statistical subtleties aside, this last result 
is obtained by expanding the transmittance in its lowest 
order in L x /£. 
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FIG. 11. Scheme of the use of the calculation of 1-d trans- 
mittance to obtain the transport properties of a 2-d system. 

A similar reasoning can be applied to higher dimen- 
sional systems (d > 2) with hyper-cubic shape. The 
conductance can be calculated with Eq.(||) by attaching 
two one-dimensional leads. The metal-insulator transi- 
tion is possible because the curve £(L) = 2£(L/a) d ~ 1 , 
valid for weak scattering, bends down to an asymptotic 
value £ (L) — > ^ . This effect manifests itself in a univer- 
sal behavior. One considers the conductance of samples 
of size L d = (Na) d with Anderson disorder with two 
one dimensional leads attached as shown in Fig. 1 1 . The 
adimensional conductance g^ = ~^G R P L (L) is calculated 
using Eq.(j^). For each dimensionality, all sizes and dis- 
orders strength scale into a single scaling curve 
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FIG. 12. Numerical scaling function j3 = {L/gr,) 
(Api/AL) evaluated from the adimensional conductance 
ql of samples of 2-d and 3-d system size L d . We use 
the Anderson model at ef = 0. Data for disorder 
W/V = 4, 5, 6, 7, 8, 10, 14 and sizes up to L 2 = 20 x 20a 2 
collapse in the 2-d curve of the filled circles. Data for disor- 
der W/V = 6, 8, 10, 14, 15.5, 16.5, 17 and 18 and sizes up to 
L A = 6 x 6 x 6a 3 collapse in the 3-d curve of the empty circles. 
In 3-d the critical conductance is g c — 3.2 and corresponds to 
a critical disorder W c /V ~ 16. 



D. The Aharonov-Bohm effect and non-local 

fa=lL*2L. (105) properties. 

9l aL 



This is shown in Fig.[l2] for dimensions d = 2 and d = 3 
as a function of ln[gi]. 

The arrows indicate the directions in which moves 
as the size of the system increases. The fact that f3d=2 
(full circles) flows from /3a = towards — oo indicates 
that, for any positive disorder all states are localized. 
In contrast, in a [3d=3 the curve can flow either towards 
02 = 1, if the disorder is weak enough (W < W c ), or to- 
wards /?2 — * — oo for disorder above certain critical value 
(W > W c with W c ~ 6.5V). This shows that the MIT 
is a critical phenomenon as any other thermodynamical 
phase transition. 

The point we want to stress here is that the "non- 
invasive voltage probes" of the Landauer picture rescue 
a "local" meaning for the conductance and emphasizes 
the non-extensive behavior introduced by quantum inter- 
ferences. The presence of decoherent processes or volt- 
meters on a scale "break" the quantum conductance 

into an incoherent sum of jf- separate pieces where the 
transport is coherent and hence described by Eq. (||). 
The total conductance becomes 0>r^l{.L x ) = ^G i R p L (L ( f ) ) 
recovering the extensive Ohmic behavior. 



Setting the Zeeman effect aside, the main consequence 
of a magnetic field is to affect the wave function's phase 
in 




In the Peierls substitution, appropriate for discrete mod- 
els, this is achieved by modifying the coupling V <— > 
Vexp[i4>]. 

When we have a systemtll as that shown in FigJT^, 
the application of the magnetic field will result in an os- 
cillation of the conduct anceESl due to the interference of 
the probability amplitudes propagated through different 
branches of the system. 




FIG. 13. Two tight binding representations of path loops 
affecting the conductance in a non-local fashion as discussed 
in the text. 
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Let's assume that in absence of field these pathways su- 
perpose in the form of maximally destructive interference 
in the forward direction and have constructive probabil- 
ity of return Pq,o (this is localization) 



model where the time is parametrically associated with 
the width represented in Fig.|l5|. 



■Po.o 



i(paths) 



U Q U J Q 



(106) 



G(2e 2 /h) 



The first term in the right side is the contribution of 
uncorrelated pathways. The second term gives the in- 
terference effects, it is affected by the magnetic field and 
contains the non-local effects intrinsic to the Schrodinger 
equation. In presence of a gauge field A, destructive in- 
terference is modified producing an increase of the tpms- 
mittance. This is a positive magnetoconductanceLJ, a 
characteristic effect of coherent systems at very low fields. 
If the inelastic events along the pathways start contribut- 
ing to the loss of phase memory, the field effect becomes 
progressively weaker until it eventually disappears. As 
we will see bellow, a very strong magnetic field causes 
interferences in the short distances (toward the Landau 
levels) therefore giving a new tendency toward localiza- 
tion. 




FIG. 14. Quantized conductance as a function of channel 
width. 




E. Conductance quantization in 2-d 

Extending what is observed in a 1-d ordered chain, 
where transmittances are either or 1 depending on the 
energy of the injected excitation; when one considers the 
l + -d (strips) and real 2-d systems, we have to account for 
the different propagating modes. As the available energy 
increases, new lateral modes becomes available increasing 
the total current. According to Eq. (^7|) the conductance 
is bounded by 

e 2 

Gr,l <-mm[ L M (£F) , R M (EF) }. (107) 

For perfect transmitting samples, a situation that 
requires order and reflection symmetry ( L M( £F ) = 
R M( ep )), Tjj is either 1 or and one obtains the con- 
ductance quantization in integer multiples of e 2 /h. This 
effect is spectacularly verified in specially tailored nano- 
structures. The conductance's discrete steps are shown 
in FigQ 

Such an effect is restricted to artificial devicesc3. The 
same effect occurs when there is a metallic bridge be- 
tween two electrodes and these are pulled apart. The 
bridge stretches becoming progressively thinner until its 
breakdown. This process is monitored by the carriers, 
which, as the bridge becomes narrower, loose propagat- 
ing channels below the Fermi energy. As a consequence 
the conductance shows clear steps such as those shown 
in in Fig.[l4| The adiabatic approximation, justified by 
the high electronic speeds as compared to any modifica- 
tion of the contact structure, allows the use of a simple 



FIG. 15. Geometry of neck formed when a nanowire is 
pulled lengthwise. As the transverse dimension shrinks con- 
ductance channels are excluded as shown in the lower panel 
where we show a scheme of the transversal modes energies as 
a function of position. 



F. Hall Effect 

In the classic theory of the Hall effect the external mag- 
netic field produces a coefficient Rh = —1/ne which is 
negative if the carriers are electrons and is inversely pro- 
portional to the electron density. In a two dimensional 
system, we can evaluate the Hall effect at low fields by 
using a simple model that has two l-d contacts to in- 
ject and measure the current and two additional ones to 
measure the Hall voltage. These are shown in Fig.O 




FIG. 16. Geometrical representation of the model used to 
evaluate the Hall effect. 
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Since the chemical potentials and voltages are related ^-^ 
by V p = (Xp/e, Eq.(§) becomes 'p = G Pi9 (V p - V g ). 

4 
9=1 

TTrl „ ,, , . . , . , . ,, . We can always choose one of the voltages, e.g. one of the 

When all the chemical potentials are equal, there is no , , , £ ■ ■ • \T n mu r 

c ctady^G — Y^G f o h'ch current leads as a reference setting V4 = 0. iherefore, 

^ p ' 9 4^ 9,p the problem simplihes to three linear equations 



li \ /G 2 ,i + G 3 ,i + G 4 ,i -Gi, 2 -Gi, 3 W Vi \ 

b I = I — G2,l G12 + G32 + G4 ; 2 — G2,3 II V2 

I3 / \ — G3,i — G3 : 2 Gi^ + G 2: 3 + G 4j 3 J y V3 y 

In short 1^ = y\ -M^V 'j, where we invert the matrix of conductances M obtaining 

Vi\ /Ri,i Ri, a Ri, 3 \ / li \ 

V2 I = I R2,l R2,2 R2,3 II b ■ 
V3 / \ R3,l R3,2 R3,3 / \ b / 

Let's see one of its elements, 

R _ (G 2 ,3 + G2,4 + G 2 j)(G 3i i + G3 i2 + G34) — G 2i 3G3 ; 2 

1,1 ~ det[M] ' 



The voltage between 2 and 3 is related to the total cur- 
rent 

V V 2 - V 3 

r<H — 7- — j — r< 2 ,i — K3.1 

I ll 

which will be zero in absence of a magnetic field. Oth- 
erwise, the dependence on magnetic field implicit in the 
quantum mechanical evaluation of R2.1 and R3.1, gives 
the Hall resistance. 

There are certain symmetry relations that should be 
obeyed on the basis of the so called Onsager-Casimir re- 
lations; in the linear response regime the Hall resistance 
obeys 

R H (B) = -R H (-B), (108) 

where B is the applied field. The previous relation in- 
volves the minus sign because the Hall voltage is reversed 
in the field. Nevertheless, if one simultaneouly exchanges 
the voltage probes there is no change. On the other hand 
the ordinary longitudinal resistance does not change un- 
der the reversal of the field. The symmetry relation can 
be summarized as 

R XX (B) = R xx (-B), (109) 

and 

R xy (B) = R yx (-B). (110) 



These symmetries are well satisfied for an homogeneous 
sample. On the other hand, for inhomogeneous samples 
no specific symmetry is experimentally found upon re- 
versing the field. This is due to irregular current patterns 
occuring in the sampple mixing both xx and xy compo- 
nents the first one being symmetric, and the second an- 
tisymmetric, in the field. A more general relatioship for 
a four probe device can be given as follows: if one de- 
fines the resistance R m n,ki as that taken when measuring 
the voltage between probes k, I while putting current / 
between probes to, n (imput at to output at n) 

r Vk ~ V| nm 



The reciprocity relations are then expressed as 

RfcZ, mn(-B) = Rmn,fci( — B) . (H2) 

The Integer Quantum Hall Effect IQHE0 follows from 
the use of Eq. ^ in the condition that strong mag- 
netic fields confine the propagating electrons to the sam- 
ple boundaries. The appearance of these new spa- 
tial selectioiu rules precludes backward scattering and 
decoherenceea, yielding the unitary transmittances re- 
sponsible for the IQHE. 
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VI. DECOHERENCE IN QUANTUM 
TRANSPORT 

A. Phenomenology 

It is obvious that because of the experimental limits 
for the coherent description introduced above one should 
observai3 important departures from those predictions. 
A first alternative to include decoherence in quantum 
transport was inspired by the Landauer's formulation. 
There, the leads, while accepting a quantum description 
of their spectra and their ability to propagate excitations, 
are the ultimate source of irreversibility and decoherence: 
Electrons leaving the leads toward the sample are com- 



pletely incoherent with the electrons coming from the 
other leads (see Fig. [l?]). In fact, it is obvious that a 
wire connected to a voltmeter, by "measuring" the num- 
ber of electrons in it, must produce some form of collapse 
of the wave function leading to decoherence. 

Besides, no net current flows toward a voltmeter. 
Leads are then a natural source of decoherence which can 
be readily described in the Landauer's picture if one uses 
the Landauer conductances together with the Kirchhoff 
balance equations. This fact was firstly realized by M. 
Biittikcrcil. The procedure, resulting from the applica- 
tion of Kirchhoff law to each lead, is called the Landauer- 
Biittiker equationta. Let us see how it works for the case 
of a single voltmeter. In matrix form 





-( 


• 





G<£,£ 

G_R,L 



Gl,0 
Gfl,^ 



Gl,r 
[G^,_r + Gl,r) 




(113) 



r 



Here, the unknowns are II, \r and = Sfi^/e. The sec- 
ond equation must be solved with the voltmeter condition 

U = o 

6 C 

= -T^ L (Sfj,^ - 5{i L ) + -Tr^(<W - Sfj, a ), (114) 
yielding us Sfi^ to be introduced in the third equation 

C 6 

l-R = ~^Tr,l {&fi>L - <Vr) + -T Rt< f,(Sfj, L - 5n<j,) (115) 



to obtain the current at the voltage source 



h = t t R,l(S^l 
h 



S^r) 



with 



Tr,l — Trx + 



(116) 



(117) 



The first term can be identified with the coherent part, 
while the second is the incoherent or sequential, i.e. the 
contribution to the current originated from particles com- 
ing from the voltmeter. This corresponds to an effective 
conductance of 



3R,L 



1R,L 



(G 



R,<p 



G<A,L ) J 



(118) 



which can be identified with the electrical circuit of Fig. 
|l8| . This classical view clarifies the competition between 
coherent and incoherent transport. What this circuit 
does not hint at is that in Quantum Mechanics one can- 
not modify one of the resistances without deeply altering 
the others. 



reservoir 




sample 










lead 





reservoir 



reservoir 



FIG. 17. Representation of a three probe measurement. 
The volt-meter may be strongly coupled and is a source of 
decoherence. 




R 



FIG. 18. Classical circuit representation of the 
non-classical system with quantum coherent and incoherent 
transport. The coherent component is a direct resistance be- 
tween left (IS) and right (R) electrodes. Quantum mechanics 
makes these effective resistances interdependent. 
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So far with the phenomenology. The next important 
step is to connect these quantities with actual model 
Hamiltonians. This connection was mad&,explicit by the 
contribution of D'Amato and PastawskiE3(DP). 



B. The D'Amato-Pastawski model for decoherence. 



The DP model refers to a simple way to account for the 
infinite degrees of freedom of the thermal bath or the elec- 
tron reservoirs. This follows from our general approach: 
to obtain exact solutions to simple problems instead of 
finding approximate solutions of complex problems. Let 
us first review the basic mathematical background that 
made possible the selection of a simple Hamiltonian that 
best represents the complex sample-environment system. 
The objective was to use its exact solution in the Lan- 
dauer's transport equation. Here, we describe the DP 
model for decoherence and show how it applies to a sim- 
ple resonant tunneling system. Consider the sample's 
Hamiltonian 



JV 



JV 



n = J2{Ei\i){i\+ £ V^[|*)01 + |j>(i| 

J=l(#i) 



1=1 



(119) 



where i and j indicate sites on a lattice. Notice that in- 
teractions are not restricted to nearest neighbors. How- 
ever, for the usual short range interactions, the Hamilto- 
nian matrix has the advantage of being sparse. The local 
dephasing field is represented by 



(120) 



where *T = %/{2t^) and we consider for simplicity only 
two one-dimensional current leads L and R connected at 
sites 1 and N respectively 



lea.cU 



± R = -i I 



T II) (II 



R 



r \n) (iV|) 



(121) 



We see that the 1 st site has escape contributions towards 
both, the current lead at the left, L Ti, and the inelastic 
channel associated to this site, The on-site chemical 
potential will ensure that no net current flows through 
the latter channel. 



C. The solution for incoherent transport 

The adimensional conductances are just the transmis- 
sion probabilities. In principle, they can be computed 
in terms of the Green's function according to Eq. (p^). 
To simplify the notation we define the total transmission 
from each site as 



JV+i 
= ]T T jti 

3=0 



AttN 1 L ri 



for i = 
for 1 < i < N 
\ for i = N + 1 



(122) 



The last equality follows from the optical theorem. The 
balance equation becomes 



JV+1 



U = = - {l/gi) Sfxt + £ V , 



(123) 



3=0 



where the sum adds all the electrons that emerge from a 
last collision at other sites (j's) and propagate coherently 
to site i. These include the electrons coming from the cur- 
rent source i.e. Ti^SfiL and the current drain. However, 
since we refer all voltages to the last one, Ti^SfiR, = 0. 
We remark that here we do not exclude the i-th site from 
the summations in Eqs. (123) and (|12z 



While this has 
no consequences in the steady state, they become rele- 
vant in the time dependent formulation. To fix the phys- 
ical interpretation we emphasize that the last term ac- 
counts for the electrons that emerging from a dephasing 
collision at site j will propagate coherently to site i where 
they have a dephasing collision. The first term accounts 
for all the electrons that emerge from this collision on 
site i to have a further dephasing collision either in the 
sample or in the leads. The net current is identically 
zero at any dephasing channel ( "lead" ) . The other two 
equations are 

JV 

\ L = -\ = - (l/ 9i ) Sfx L +Y,T itj 5n3, (124) 

j=o 

JV+l 

\ R = I = - (l/gi) 6^ + £ TijSfij. 

3=0 

Here, we need the local chemical potentials. They can be 
obtained from Eq. (123). In a compact notation, these 
coefficients can be arranged in a matrix form which ex- 
cludes the leads that are current source and sink 



W 



T 2,l 



3,1 



^1,2 ^1,3 ••• 7\,JV 

^2,2 — l/ff2 T 2j 3 ••• T 2: AT 

^3,2 r 3)3 — l/g 3 ■■■ T 3i jv 



(125) 



Tn.: 



T N . 



JV 
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from which the chemical potential in each site can be 
calculated as 



5 f i i = '£[W- 1 ]..T j , 8» 



(126) 



Replacing these chemical potentials back in Eq. ( |123| ) 
the effective transmission can be calculated 

N N 

j=l i=i V Ll > ) 



coherent 



incoherent 



The right side contains two contributions: the first 
comes from electrons that propagate quantum coherently 
through the sample, the second contains the incoherent 
contributions due to electrons that suffer their first colli- 
sion at site i and their last at site j. 

Until now the procedure has been completely general, 
there is no assumption involving the dimensionality or ge- 
ometry of the sample. The system of Fig.[l9] was adopted 
by DP only because it has a simple analytical solution 
for various situations ranging from tunneling to ballistic 
transport. We summarize the procedure for the linear 
response calculation. First, we calculated the complete 
Green's function in a tight binding model. With the 
Green's functions we then evaluate the transmittances 
between every pair of sites in the sample (i.e. nodes in 
the discrete equation) and write the transmittance ma- 
trix W. Then, we solve for the current conservation equa- 
tions that involves the inversion of W. 




r 



Hi fll fll ft 



FIG. 19. Pictorial representation of the 
D'Amato-Pastawski model for the case of a linear chain. 

What are the limitations of this model? A conceptual 
one is the momentum demolition produced by the local- 
ized scattering model. Therefore, we go directly from 
quantum ballistic description to a classical diffusive one. 
To describe the transition from quantum ballistic to clas- 
sical ballistic, one should modify the model to have the 
scattering defined in momentum or energy basis. 

The other aspect is merely computational. Since the 
resulting matrix W is no longer sparse, this inversion is 
done at the full computational cost. A physically appeal- 
ing alternative to matrix inversion was proposed in DP. 



The idea was to expand the inverse matrix in series on 
the dephasing collisions, resulting in 

i i j 

+ /] I H..U.I Li. + ■■■ (128) 

i j I 

The formal equivalence with the self-energy expansion of 
Matsubara and Toyosawa in terms of locatorsc3 or local 
Green's function justifies identifying gi as a loeator for 
the classical Markovian equation for the density^ gener- 
ated by the transition probabilities. 

As an example of use of Eq. (12S), one can readil y ap - 
ply it for the model of Section VIA obtaining Eq.(117). 
This also constitutes the basis for a perturbative method 
of calculating the conductance with a substantially re- 
duced computational cost. This strategy has been ap- 
plied with considerable success to explain the stability 
of the Quantum Hall Effecto against scattering and de- 
phasing. 



Notice that Eq. (128) can also be rearranged as 



N 



T, 



R,L 



Tr.l + J2 Th >* x 



T, 



total = coherent + totalx 1 st collisionx coherent 

(129) 

The summation on the right hand side has the structure 
of the Dyson equation. It is a Bethe-Salpeter equation. 
This is graphically represented in Fig.^. We notice that 
according to the optical theorem gi — 2?:hT^Ni, while 
both transmittances entering the vertices in the figure 
are proportional to l/r^ the whole product involved in 
the vertex is proportional to the dephasing rate. The ar- 
rows make explicit that transmittances are the product 
of a retarded (electron) and an advanced (hole) Green's 
function. To describe diffusive scattering by impurities 
one must replace the whole electron-hole pair by its en- 
semble average calculated in the ladder approximatiorL-nr 
eventually include their quantum corrections (see Ref.Cj) . 

Since the transmittance is essentially a two-particle 
Green's function, now it is a particular case of the Bethe- 
Salpeter equation which is solved in the "ladder" ap- 
proximation whichpEpr the present model is exact (see 
appendix B in Rcf.c3 for a detailed formalization of this 
picture with the Keldysh formalism). In summary, de- 
pending on the approximation used for the density prop- 
agator, the self-consistent solution contains the metal- 
lic transport, the weak-localization corrections and even 
the thermally activated hopping regime. It also, can be 
viewed as a self-consistent Born approximationcZI for the 
electron-phonon interaction. One might also want to in- 
troduce other forms of interaction, such as those which 
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conserve momentumi 
point here. 



but we will not extend on this 




gi = Xi/dnhNi ) 



FIG. 20. Feynman Diagram for the Dyson Equation of the 
transmittance. It is equivalent to a particle-hole Green's func- 
tion in the ladder approximation where the usual interaction 
rung connecting the two propagators is here represented by a 
dot. 

Many of the results contained in the D'Amato- 
Pastawski paper for ordered and disordered systems were 
later extended in great detail in a series of papers by S. 
Datta and are described in his didactic bookE3. In the 
next section, we will illustrate how the previous ideas 
work by considering again our reference toy model for 
resonant tunneling. 



D. Effects of decoherence in resonant tunneling 

If jsie consider the "sample" to consist of a single 
stateO. If we choose to absorb the energy shifts into 
the site energies E = E + A (e) , the Green's function 
is trivial 



G£ (e) 



1 



e-E + i( L T + R T 



(130) 



By taking the T's independent on e in the range of inter- 
est, we get the "broad-band" limit. We drop unneeded 
indices and arguments for the time being. From this 
Green's function all the transmission coefficients can be 
evaluated at the Fermi energy 



(131) 



From the energy dependent transmittance we obtain the 
total transmittance 
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l t + r r 



(132) 



The first term in the curly bracket is the coherent contri- 
bution while the second is the incoherent one. We notice 
that the effect of the decoherence processes is to lower the 
value of the resonance from its original one in a factor 



( L r+ R r) 



*r)' 



(133) 



In compensation, transmission at the resonance tails be- 
comes increased. 

It is interesting to note that if the resonant level lies be- 
tween fi a + eV and fi D and provided that the voltage drop 
eV is greater than the resonance width ( L T + R T + ^T), 
we take Tr^e, eV) ~ Trl(e) and we can easily com- 
pute the non-linear response. Notably, one gets that the 
total current does not change as compared with that in 
absence of decoherent processes, i.e. 

lf=/ T RtL {e)de 



fJ-o+eV 



T° R , L {e)de 



= 47T R T ~7y~ — — TT^TT" L r. 



( L r+ «r) 



(134) 



Thus, in this extreme quantum regime, the decoherence 
processes do not affect the overall transport. 

In spite of the complexity of the general problem of de- 
coherence in mesoscopic systems we extract an important 
lesson from the case of resonant tunneling with decoher- 
ence solved above within the DP model. The inclusion 
of external degrees of freedom has three effects 

1) It broadens the resonance relaxing energy conserva- 
tion. 

2) The integrated intensity of the elastic (coherent) 
peak is decreased. 

3) The inelastic contribution came out to compensate 
this loss and maintains the value of the total transmit- 
tance integrated over energy. 

4) Taking this system as representative of those whose 
spectra are strongly quantized, i.e. with well defined, iso- 
lated resonances, one may state that these systems are 
quite stable against decoherence. 



VII. THE SOLUTION OF TIME DEPENDENCE 

First we remember that the Green's function gives the 
possibility to write non trivial initial values in the space- 
time coordinate Xj = (jj,tj) described by Source (Xj) ■ 



Therefore, any wave function injected at space time 
Xj = (r j , t j ) by an arbitrary source is propagated by 
the retarded Green's function 



<p(X 2 ) = ih / G R (X 2 ;X J ) Vsourco {X J )dX J 



(135) 
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Here, instead of using a discrete spatial index, we 
adopted the notation G^ 2 r .(t 2 -tj) — > G R (X 2 ; Xj) and 
J2 r f dij — > f dX,. To inject particles with some defini- 
tion in momentum and/or energy we need to introduce 
precise correlation between space and time variables. We 
resort to the complementary equation 

<p*{X t ) = -ih [ ^ urce (X fe )G A (X fe ;Xi)dX k , (136) 



which we can use to find how the density matrix of par- 
ticles depends on both correlated initial conditions 

[<p*(X 1 )<p(X a )] = 7? J J G R (X 2 X 3 ) 

x [cp SOUI UX J )cp: omcc (X k )}G A (X k X 1 )dX i dX k (137) 

The usual density function is defined by taking X\ = 
X 2 i.e. p{Xi) = tp*(X 1 )(p(X 1 ). This is nothing 
else but the Schrodinger equation written in a gen- 
eral form that allows arbitrary boundary conditions 
[ipsourcc(Xj)(pl OUICC (Xk)} . One is used to the initial 
conditions at space-time coordinate Xi of the form 
Psource (Xi) = [if source (Xi)<p* OUIce (Xi)] which allows pre- 
cision in the specification of position and time at the price 
of absolute uncertainty in momentum and energy. The 
formalism allows one to program an uncertainty trade- 
off to approach better to a semiclassical initial condition 
with energy and momentum known up to some precision. 

The basic idea to get the physics of time dependent 
phenomena from Eqs. (pL3 7 | ) n Qy from its generalization 
in Quantum Fields TheoryE3Ea is to recognize that in 
any Green's function G(tj;tk), a macroscopically ob- 
servable time is t = \ [tj + tt] whose Fourier transform 
is an observable frequency lo. Meanwhile, system's mi- 
croscopic energies, e, are Fourier transforms of internal 
time differences tj — %. By operating carefully with 



this concept in Eq. (137), the dynamical transmittances 
T(r/,rjj£, w) = Tjs(e, cu) with n a position in channel i 
and Yf a position in channel / are obtained. The manip- 
ulatiorL_Qf the time integrals is quite subtle and we refer 
to Ref.t2l for details. The basic result is 

T f 4e,u,) = 2T f Gf l (e+hu;)2r i Gf J (e~hu;) 7 (138) 

and is consistent with the steady-state Fisher-Lee equa- 
tion. From this new formula we can evaluate the proba- 
bility that the wave packet of mean energy e propagates 
from i to / in a time t as 



Tf,i(e,t) = J T A4 (e,cj)cxp[-icjt] 



(139) 



This time dependent transmittance was incorporated in 
a generalization of the Landauer formula for time depen- 
dent phenomena. 

Since the spectrum is continuous we can keep the low- 
est order in the frequency expansion and obtain 



1 — iwrp 



(140) 



Clearly, such approximation would give an exponential 
for the propagation dynamics which does not describe 
the rKery short time regime. According to Eq.(4.5) in 
Ref.Ej, the typical propagation time Tp between points 
Xi and r^- results 



ih 

Tp = 

2 



GfA^Gf^-' + GtAe^Gfjie)- 1 



(141) 



de 

ihd_ Gf/e) 
'2ds n G^(e) 



This was evaluated in various simple systems in Ref.tii 
There, for ballistic metals with velocity v e one gets 
Tp = \y f — Yi\ / v e . For diffusive metals where impurity 
collision at a rate 1 /r determines the diffusion constant 
D e = v 2 £ t /2 it results t p = |iy - r,| 2 /(2D e ). 

As a striking example, we mention the "simple" case 
of tunneling through a barriered of length L and height 
U exceeding the kinetic energy e of the particle. One 
gets 



t p = 



(142) 



which, within our non-relativisticjiescription, can be ex- 
tremely short (even super lumina£3!) provided that the 
barrier is high enough. 

In a double barrier system, in the regime of resonant 
tunneling, the above contribution becomes negligible and 
the propagation time is then fully determined by the life- 
time inside the well. In fact, by using the functions of 
the previous Section, one gets a delay of 



T P = 



2( L r+ «r) 



(143) 



which readily limits the admittance associated to the 
device to G(w) = G(0)/(1 — iuiTp). This is in fan- 
agreement with the experimental resultsH Tunneling 
times can also be calculated in more complex situa- 
tions such as disorderedEilsystems and situations with co- 
herent capacitive effects and phonon interactionsE3 and 
decoherentE^processes. . 



VIII. FINAL REMARKS. 

We hope these lectures have been able to convey, at 
least partially, the essence of our message: In our at- 
tempt to understand the quantum world it is possible, 
and convenient, to make very simple models of nature. 
Those models, though not necessarily complete, can show 
us many effects that surprise our classical intuition and 
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which would have been obscured by other more "com- 
plete" descriptions. Every one of these surprises can give 
a new twist in the experimental research, which is indeed 
the ultimate truth, and eventually become the source of 
an innovative application. 

In our toolbox to extract information from the pro- 
posed models the most important elements have been the 
decimation method and the Green's functions technique. 
Both are intimately connected with the renormalization 
group concept. As such they arc fundamental in generat- 
ing the models themselves, as they help us to identify the 
relevant variables. Simultaneously, our understanding of 
the quantum world has benefited from fresh approaches 
to long standing problems as that of the tunneling time. 
In the process of computing transport properties we have 
also learnt the deepest principles of statistical mechanics 
as to what is the origin of macroscopic irreversibility and 
established the conditions for the validity of the ergodic 
hypothesis. The fundamentals of decoherent processes in 
transport and their consequences are currently of intense 
interest. The lines of thought we have sketched here, 
have turned out to be very productive in this direction, 
and we hope they will continue being fruitful through the 
action of our students. 
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